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Abstract. We observed the main core F of the p Ophiuchi cloud, an active star-forming region located at d ~ 140 
pc, using XMM-Newton with an exposure of 33 ks. We detect 87 X-ray sources within the 30' diameter field-of- 
view of the EPIC imaging detector array. We cross-correlate the positions of XMM-Newton X-ray sources with 
previous X-ray and infrared (IR) catalogs: 25 previously unknown X-ray sources are found from our observation; 
43 X-ray sources are detected by both XMM-Newton and Chandra; 68 XMM-Newton X-ray sources have 2MASS 
near-IR counterparts. We show that XMM-Newton and Chandra have comparable sensitivity for point source 
detection when the exposure time is set to ~ 30 ks for both. We detect X-ray emission from 7 Class I sources, 26 
Class II sources, and 17 Class III sources. The X-ray detection rate of Class I sources is very high (64 %), which 
is consistent with previous Chandra observations in this area. We propose that 15 X-ray sources are new class 
III candidates, which doubles the number of known Class III sources, and helps to complete the census of YSOs 
in this area. We also detect X-ray emission from two young bona fide brown dwarfs, GY310 and GY141, out of 
three known in the field of view. GY141 appears brighter by nearly two orders of magnitude than in the Chandra 
observation. We extract X-ray light curves and spectra from these YSOs, and find some of them showed weak X- 
ray flares. We observed an X-ray flare from the bona fide brown dwarf GY310. We find as in the previous Chandra 
observation of this region that Class I sources tend to have higher temperatures and heavier X-ray absorptions 
than Class II and III sources. 
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C3 ■ 1. Introduction al. 1994; Kamata et al. 1997; Tsuboi et al. 2000), and 

Chandra (Imanishi et al. 2001a; Imanishi et al. 2001b; 

"O ' Y ° Ung l0W " maSS StarS are kn0wn t0 be ubic l uitous X " m y Imanishi et al. 2003). To investigate the YSO X-ray emis- 

■ emitters (Montmerle et al. 1993; Feigelson and Montmerle sion in thc dcnsc C0 core F of the p Ophiuchi cloud (Loren 

g ! 1999). Their X-ray luminosities are ~ 10 4 - 10 5 times et aL 1990 ) ) we took with XMM-Newton a 33 ks exposure 

higher than that of the present-day Sun. YSOs are classi- of this reg i n, as part of the EPIC Guarantee Time pro- 

fied from IR and sub-millimeter spectral energy distribu- gramj and we report here the results of this observation, 

tions (Andre & Montmerle 1994) into four classes: Class In Sect we pre sent observation and data analysis, in- 

and Class I protostars, Class II (classical T Tauri stars) eluding source detection, extraction of X-ray light curves 

and Class III (weak-lined T Tauri stars) sources, inter- and X-ray spectra from detected sources. We compare our 

preted as a time evolution sequence from Class (age XMM-Newton results with previously obtained Chandra 

~ 10 4 yrs) to Class III (age £ 10 6 yrs). It is generally ac- results in this region in Sect W e discuss the IR prop- 

cepted that enhanced magnetic activity of YSOs produces erties of detected sources in Sect. ED We present the X-ray 

X-rays emitted by high temperature plasmas as a result of detections of young bona fide brown dwarfs in Sect. ® and 

heating by magnetic reconnection events, which is called the X - ra y properties of Class I, II and III sources in Sect, 

the solar paradigm because of the analogy with the X-ray rg] Fma n Vj we summarize all results in Sect. [7| 
emission of the Sun (e.g. as seen with Yohkoh). 

The p Ophiuchi cloud is a well-known nearby ac- 
tive star-forming region at d ~ 140 pc, which has been 



2. Observation and data analysis 



observed over twenty years by practically all X-ray ob- The main (densest) core F of the p Ophiuchi cloud 

servatories: Einstein (Montmerle et al. 1983), ROSAT was observed with XMM-Newton for 33 ks on February 

(Casanova 1994; Casanova et al. 1995; Grosso et al. 1997; 19, 2001 pointing at a J20 oo = 16 h 27 m 26.0 s and (5 J2 ooo 

Grosso et al. 2000; Grosso 2001), ASCA (Koyama et = -24°40'48.0". XMM-Newton is the X-ray astronomi- 
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cal observatory of the European Space Agency which was 
launched on December 10, 1999. XMM-Newton has three 
co-aligned X-ray telescopes with 6"-FWHM angular reso- 
lution at 1.5 keV and with large total effective area up to 
5000 cm" 2 at 1 keV (Jansen et al. 2001). The European 
Photon Imaging Camera (EPIC), i.e. two MOS (namely, 
MOS1 and MOS2) and one PN X-ray CCD arrays, are 
placed on each focal plane of the three X-ray telescopes. 
EPIC provides imaging capability over a wide field-of-view 
of 30' diameter, and spectroscopic capability with moder- 
ate spectral resolution (~ 60 eV at 1 keV) in the 0.2—10.0 
keV band. In this observation, EPIC was operated with 
the medium optical blocking filter and in the Full-Frame 
mode. 

We reduced the data using the XMM-Newton Science 
Analysis Software (SAS version 5.4). The photon event 
lists were obtained by running emchain and epchain. We 
rejected bad pixels, and selected single, double, triple, and 
quadruple pixel events for MOS, and single and double 
pixel events for PN using evselect. We excluded the time 
intervals where the X-ray counts from whole field of view 
are extremely high, above 4.4 counts s" 1 for MOS1 and 
MOS2, and 20 counts s _1 for PN, indicating a high level of 
irradiation by solar protons, and hence a high background 
level on EPIC. The sum of the final good-time-intervals 
are 31 ks for MOS1, 32 ks for MOS2, and 29 ks for PN. 

2.1. Source detection 

Fig-fflshows the resulting EPIC image where the MOS1, 
MOS2, and PN images were corrected for vignetting and 
co-added. Red, green, and blue colors code the X-ray 
events in the 0.5-1.0 keV, 1.0-2.4 keV, 2.4-8.0 keV en- 
ergy bands, respectively. The image is smoothed to en- 
hance the contrast. 

We perform the source detection on the individual de- 
tector X-ray images, using the standard source detection 
method of SAS. This method includes 4 steps, heritage of 
the Einstein and ROSAT source detection algorithms. 

1) The source candidates are obtained by the sliding 
cell method in eboxdetect. A small size cell surrounded 
by a larger size cell is moved over the X-ray image. The 
background counts are estimated from the area surround- 
ing the small size cell. If the counts in the small size cell are 
significantly larger than in the background, the cell is con- 
sidered to have excess counts. After sliding the cell over 
the whole region, we obtain "islands" of excess counts, 
which are considered as source candidates. 

2) The local background image is created by 
esplinemap. The source candidates obtained in step 1 are 
subtracted from the X-ray image, and a spline fit is applied 
to this "swiss-cheese" image to obtain the background im- 
age. 

3) The sliding cell method is again performed with 
eboxdetect to obtain more reliable source candidates. 
This time, only the small size cell is moved over the X-ray 



image. The local background counts are estimated locally 
using the background image obtained in the previous step. 

4) A two-dimensional fit is performed on the X-ray 
source candidates in emldetect using a model where the 
X-ray counts come from a point source, contamination 
by nearby point sources, and the background. The point 
spread function is used as the model function of the point 
source. The local background created at step 2 is used as 
the background model. Through this task, the significance 
of the source candidates is checked precisely and source 
counts are obtained accurately. 

This method can deal with images in different energy 
bands and different detectors simultaneously. We use im- 
ages in the 0.3—2.0 keV and 2.0—8.0 keV energy bands, 
and we deal with the MOS1, MOS2, and PN X-ray images 
simultaneously to obtain good statistics. For the areas cov- 
ered by only two detectors or one detector, we can analyze 
only the images of the corresponding detector set. 

As a result of this protocol, we detect in total 87 X-ray 
sources with a likelihood threshold above 12, correspond- 
ing to the significance of ~ 4.4 a for Gaussian statistics 1 , 
within the 30' diameter field-of-view of the EPIC de- 
tectors. The positions of the detected 87 X-ray sources 
are indicated in the lower panel of Fig. ^ Table 1 gives 
the coordinates, count rates, and the hardness ratios of 
the detected sources. For convenience, these sources are 
designated here ROXN-n (for "Rho Oph X-ray sources, 
Newton, number n") 2 . The hardness ratios are calculated 
after the formula HR = (H - S)/(H + S), where S and 
H are counts in the 0.3-2.0 keV, and in the 2.0-8.0 keV 
energy band, respectively. 3 

The detected sources are cross-correlated with: X-ray 
source catalogs from Einstein (Montmerle et al. 1983), 
ROSAT (Casanova et al. 1995; Grosso et al. 2000), ASCA 
(Kamata el al. 1997), and Chandra (Imanishi et al. 2001a); 
near-IR source catalogs (e.g. Barsony et al. 1997; 2MASS 
= Two Micron All-Sky Survey); ISOCAM mid-IR catalog 
(Bontemps et al. 2001); and an optical catalog (Monet 
1996). Table 2 shows the resulting cross-correlation: 62 
ROXN sources are identified with previously known X- 
ray sources, hence we find 25 new X-ray sources from our 
observation, and 68 ROXN sources have 2MASS near-IR 
counterparts. Table 2 also lists the IR classifications from 
Bontemps et al. (2001), or Andre & Montmerle (1994). 

As a result of these cross-identifications, we find 7 
Class I sources, 26 Class II sources, and 17 Class III 
sources. We present 15 previously unknown Class III can- 
didates in section 01 

1 There is a coding error in emldetect until SAS version 5.4 
(XMM-Newton-News #29, ll-Mar-2003), producing an overes- 
timate of the source likelihood. We corrected the source likeli- 
hood using the table provided in the XMM-iVewfon-Newsletter. 

2 The official naming convention for XMM-Newton sources 
is XMMU JHHMMSS.S+ / — DDMMSS. 

3 ROXN-1 was so bright that the CCD pixels on the source 
are strongly affected by pile-up. Hence HR computed for this 
source is not reliable. 
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Fig. 1. XMM-Newton observation of the p Ophiuchi cloud. The upper panel shows the tricolor image where EPIC 
MOS1, MOS2, and PN images are added together after vignetting correction. The image was smoothed with a 6"- 
FWHM Gaussian to increase the contrast. The red, green, and blue colors code X-ray events in the 0.5—1.0 keV, 
1.0—2.4 keV, and 2.4—8.0 keV band, respectively. The lower panel shows the same area with the positions of the 87 
X-ray sources detected by XMM-Newton labeled by ROXN numbers (Table 1). The dashed lines show the field of view 
of Chandra (Imanishi et al. 2001a, Imanishi ct al. 2003). The background image is a combination between DSS2-red 
(black & white), ISOCAM LW2 6.7 /im (blue), and LW3 14.3 /im (red) filters. 
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Table 1. The list of the X-ray sources detected by XMM-Newton. 
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(1) XMM-Newton source ID number, "Rho Oph X-ray sources, Newton".; (5-10) The hyphens indicate that the X-ray sources are outside the field 
of view of the instrument.; (6), (8), (10) The hardness ratio is calculated by the equation, H.R. — (H — S)/(H + S), where S and H are counts in 
the 0.3—2.0 keV band and in the 2.0—8.0 keV band, respectively. Note that due to pilc-up, H.R. is maybe not reliable for source 1. 



Ozawa et al.: XMM-Newton observation of the p Oph cloud 
Table 2. Identification of the XMM-Newton sources with known X -r&y sources i md optical/infrared sources. 
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9.87 


1.92 


0.86 


GY410 


III 


77 




55 


- 




F34 


- 


- 


16275996-2448193 


10.81 


9.76 


9.27 


1.05 


0.49 


BKLT162800-244819 


nlll 


78* 
















16280464-2434560 


16.92 


13.18 


10.97 


3.74 


2.21 


GY463 


nlll 


79* 
80* 
















16280478-2437100 


13.62 


11.45 


10.34 


2.17 


1.11 


GY465 


nlll 


81* 
































82* 
83* 
















16280810-2445121 


14.12 


12.86 


12.15 


1.26 


0.71 


GY478 


nlll 


84* 
































85* 
86 




































33 










16282373-2441412 


11.73 


9.83 


9.01 


1.90 


0.82 


BKLT162823-244140 


nlll 


87 






34 




F37 






16282516-2445009 


6.72 


6.60 


6.51 


0.12 


0.09 


HD148352f 





(1) XMM-Newton source ID number, "Rho Oph X-ray sources, Newton". The asterisks indicate new X-ray sources found in this observation; 
(2 — 9) The identification with source catalogs of (2) Einstein (Montmerlc et al. 1983), ROX— ; j" Grosso (2001) shows that one of these Einstein 
observation has detected an X-ray source associated with the YLW16 IR source group; (3)ROSAT PSPC (Casanova et al. 1995), ROXR1 — ; 
{A)ROSAT PSPC (Casanova 1994), ROXR2-; (5)ASCA (Katama et al. 1997), ROX A—; (6)ROSAT HRI (Grosso et al. 2000), ROXR-; (7)Chandra 
(Imanishi et al. 2001a); (8)ISO ISOCAM (Bontemps et al. 2001), (9) 2MASS. The dotted lines indicate that the XMM-Newton sources have no 
known counterpart. The hyphens indicate that the XMM-Newton sources are outside the field of view of the instruments. (10 — 14)If only the lower 
limits of the values arc available, ">" is put before the values. (15) For the source naming convention, see Andre & Montmerlc (1994) and references 
therein. The source names beginning with BKLT are from Barsony et al. (1997). 0600-2053 is from the catalog USNO-Al.O (Monet et al. 1996). % 
indicate foreground stars. (16) The classifications of Class I, II, and III are from Bontemps et al. (2001) or Andre & Montmerlc (1994). The Class 
III candidates proposed in this paper are indicated by "nlll". The bona fide brown dwarfs arc also indicated in this column by "bd". 
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2.2. X-ray light curves 

We extract background-subtracted X-ray light curves for 
the ROXN sources. X-ray counts are selected from circular 
regions around the sources where 90% of the X-ray counts 
from the sources are included at 1.5 keV. 4 In crowded 
regions, we adjust the circle radius to avoid contamina- 
tion by neighbouring sources. Background X-ray counts 
for each source are extracted from annular regions with 
a 60" inner radius and a 150" outer radius around the 
sources. If the annuli contain other X-ray sources, regions 
within a 40" radius around those sources are excluded 
from the corresponding annuli. To avoid decreasing too 
much the area of the annuli by excluding nearby X-ray 
sources, the outer radius was adjusted to keep a constant 
area for the background extraction region. 

Since YSOs are known in general to exhibit X-ray 
flares, we systematically look for flares in our sources. We 
find that nearly half of the sources (40) have enough X-ray 
counts (JV > 150, where N is the sum of MOS1, MOS2 
and PN counts) to search for X-ray flares, using the follow- 
ing method. We extract X-ray light curves with a uniform 
reference time binsize. A bin-to-bin variation was then de- 
fined to be a flare if the X-ray counts in a given time bin 
are 1.5 times higher than at least one of the three pre- 
ceding bins with a statistical significance above 4 a. We 
apply this method varying the reference time binsizcs from 
1000 s to 4000 s by steps of 1000 s, which covers the typi- 
cal rise times of X-ray flares from YSOs. To plot the light 
curves, we use an adaptive binning where the time binsize 
is tuned to keep a constant number of counts inside each 
time bin. This enables us to catch any rapid time vari- 
ability such as rising phase of X-ray flares. Depending on 
the total number of counts, the time bins are tuned so as 
to have enough counts in each bin to see the overall time 
profile properly (see tuning given in Table 3). 

Fig. |21 shows a sample of background-subtracted X- 
ray light curves. The X-ray flares seen in Fig. show a 
variety of profiles, broadly characterized by a fast increase 
and a roughly exponential decay, although some profiles 
are more symmetrical, with comparatively slow rise and 
decay phases of similar durations. Contrary to previous X- 
ray observations which had limited statistics, a "typical" 
YSO X-ray flare profile cannot, in general, be reduced to 
a simple fast rise followed by an exponential decay. Such 
a complex behavior is also seen in the large X-ray flare 
database obtained in 40 ks Chandra observations of the 
Orion Nebula Cluster (Feigelson et al. 2002a). 

In our p Oph observations, X-ray flares are detected 
from 1/4 Class I, 6/13 Class II, and 3/12 Class III sources. 
These ratios are 2-3 times lower than those obtained by 
Chandra (Imanishi et al. 2001a; 14/18 for Class I + Class 
I candidates, 12/20 for Class II+III), which is to be ex- 
pected since our XMM-Newton exposure (33 ks) is three 
times shorter than that of Chandra (100 ks). The unclas- 
sified sources did not show any strong flares. 

4 See calibration files XRT1 XPSF 0005.CCF, XRT2 XPSF 
0005.CCF, and XRT3 XPSF 0004.CCF. 



Table 3. Tuning of the adaptive binning. Counts inside a 
time bin for a given source total counts. 



Source total counts 


Counts inside a time bin 


0- 500 


15 


500 - 1000 


25 


1000 - 2000 


35 


2000 - 4000 


50 


4000 - 6000 


80 


6000 < 


100 



2.3. X-ray spectra 

We now compute background-subtracted X-ray spectra for 
sources with N > 150 cts (4 Class I, 13 Class II, 12 Class 
III sources, and 11 unclassified sources). The way to ex- 
tract X-ray counts and background counts is the same 
as for extracting the light curves. 5 Since the difference 
between the responses of MOS1 and MOS2 is not large 
because they use the same type of CCD array, we add the 
MOS1 and MOS2 spectra to improve the statistics, i.e., we 
add the ancillary response files, and average the photon 
redistribution matrixes of MOS1 and MOS2 weighted by 
the exposure time of each detector. We treat the PN spec- 
tra separately. Fig.[3]shows examples of the X-ray spectra 
obtained in our observation. 

We first fit the separate PN and MOS1+MOS2 spectra 
with a single temperature thin thermal emission model 
(MEKAL) combined with an absorption model (WABS) 6 . 
For sources which do not have enough counts to yield the 
plasma metal abundance, we fix it to 0.3 solar, a typical 
value for YSOs obtained in p Oph with ASCA (Kamata 
et al. 1997) and Chandra (Imanishi et al. 2001a), and also 
in other star-forming regions. For calculations of emission 
measures and luminosities, we adopt d — 140 pc, based 
on Hipparcos measurements (see discussion in Bontemps 
et al. 2001), and for consistency with previous ROSAT 
study (Grosso et al. 2000). 

As a result, we obtain acceptable one-temperature 
spectral fits for all the sources except for ROXs20A + 
ROXs20B, SR12A-B and IRS55 (see below). The values 
of the fitting parameters are listed in Table 4. The individ- 
ual metal abundances obtained from the fits are found to 
be distributed around 0.3 solar, which is consistent with 
most existing YSO X-ray observations. If we alternatively 
choose to fix the metal abundance to 0.3 solar for all the 
sources, we also obtain acceptable one-temperature fits 
with better constraints on other spectral parameters ex- 
cept for the Class I source EL29 (ROXN-22). EL29 did 
not show X-ray flares during our observation, whereas it 

For ROXN-1, we exclude a circular region of 17.5" radius 
centered on the source to keep the spectrum from being affected 
by pile-up. 

6 We check consistency of Nu determination between 
Chandra and XMM-Newton in Appendix B. We mention also 
uncertainties of X-ray derived Nh values especially for an effect 
of metal abundances for the absorption model. 



Ozawa et al.: XMM-Newton observation of the p Oph cloud 



7 



M0S1 +M0S2 + PN, 0.3-8 keV 
ROXN-16, RS4-4/YLW1SA, Class 



1x10 2x10 3x10 

Tiaie(s) 



MOS1 +MOS2 + PN, 0.3-8 keV 
R0XN-65, GY31^, Class II 



1x10 2x10 3x10 

Tlme(s) 



; M0S1 +M0S2 + PN, 0.3-8 keV 
TROXN-73, IRS55, Class II 



1x10 2x10 3x10 

Timc(s) 



0.30 
0.25 
0.20 
0.15 
0.10 
0.05 
0.00 



M0S1 +M0S2 + PN, 0.3-8 keV 
ROXN-43, IRS43/YLW15, Class 



2x10 
Time(s) 



0.30 


"MOS1 +M0S2 + PN, 0.3-8 keV 


0.14 ; 


0.25 


^R0XN-30, WL20, Class II 


0.12 ; 


0.20 






0.10 ; 


0.15 






2 0.08 : 
S 0.06 ; 


0.10 


: i ,/ 


' lfl W : 


0.04 ■ 


0.05 






0.02 


o.on 






0.00 ' 


4 4 4 
1X10 2x10 3x10 






7M0S1 +M0S2 + PN, 0.3-8 keV 
lR0XN- - 8, WL17, Class I 



Time(s) 



2x10 
Time(s) 



MOS1 +MOS2 + PN, 0.3-8 keV 
ROXN-45, VSSC25, Class II 



K10 2x10 
Time(s) 



-MOS1 +MOS2, 0.3-8 koV 
ROXN-50, IRS47, Class I 



1x1(T 2x10^ 
Time(s) 



1 .0 
. 0.8 - 

§ 0.6 - 

o 

U 

0.4 - 
0.2 - 
0.0 L_ 



W0S1 +M0S2 + PN, 0.3-8 keV 
.ROXN-36, SR12A-B, Class II 



2x10^ 
Time(s) 



Fig. 2. A sample of X-ray background subtracted light curves of YSOs obtained with XMM-Newton showing variabil- 
ities. High background time intervals were suppressed (holes in the light curve). We use an adaptive binning to keep 
a constant count number in each bin (see Table 3). X-ray light curves of bona fide brown dwarfs, GY310 and GY141, 
are shown in Fig. [5] 



showed X-ray flares during the ASCA and Chandra obser- 
vations (Kamata et al. 1997, Imanishi et al. 2001a). We 
find 1.0 (0.80-1.3) for the metal abundance of EL29 rel- 
ative to solar, which is far from the value, ~0.3, observed 
for the quiescent level by ASCA and Chandra (Kamata et 
al. 1997, Imanishi et al. 2001a). Actually, the Ka line from 
He-like iron ions at 6.7 keV is clearly seen in the spectra 
of EL29 (see Fig- - We note a residual in the spectral 
fitting around 6.4 keV in both PN and MOS1+MOS2, 
however adding a gaussian line at this energy does not 
change the metallicity at all. This suggests that the metal 
abundances of the X-ray emitting plasma in YSOs could 
be variable. For example, such metalicity enhancement 
could be explained by photospheric evaporation produced 
by flare (Gudel et al. 2001). Then we decide to list the val- 
ues of the metal abundances obtained from X-ray spectral 
fitting in our observation in Table 4. 

Since it is often found that X-ray spectra with low 
statistics can be modeled by either a single temperature 



or a two-temperature MEKAL + absorption model, we 
also test a two-temperature MEKAL + absorption model 
for all the ROXN sources. The metal abundance is then 
fixed to 0.3 solar except for EL29. We thus obtain ac- 
ceptable two-temperature fits for all the sources, but the 
emission measures of the soft components are either ex- 
tremely large or very low for most sources. This means 
that the large absorption of the sources does not allow us 
to constrain the parameters of the soft component, so that 
one-temperature fits are sufficient in practice. However, for 
the bright Class III sources ROXs20A and ROXs20B (i.e., 
ROXN-27 and ROXN-28, which are spatially resolved but 
with too few counts, hence we use a circular area encom- 
passing both sources to obtain a spectrum), SR12A— B 
(ROXN-36) and IRS55 (ROXN-73), two-, three-, and two- 
temperature MEKAL + absorption models, respectively, 
are needed to obtain acceptable fits (see Fig.[j|J). 
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Table 4. Spectral parameters obtained from fit 



Class I 



ROXN 


NAME 




kT 


AT H 


abundance 


E.M. 


Lx 


x 2 /d.o.f 










(10 22 cm" 2 ) 




(10 52 cm" 3 ) 


(10 29 erg s- 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


22 


EL 29 




4.3(3.6-5.1) 


4.9(4.4-5.3) 


1.0(0.80-1.3) 


16.0(13.6-18.9) 


27.5 


90.0/ 63.0 


43 


IRS43/YLW15 




3.9(3.0-5.3) 


3.3(2.9-3.8) 


[0.30] 


13.1(10.2-17.3) 


17.3 


36.5/ 35.0 






f 


3.2(2.7-3.9) 


4.0(3.5-4.5) 


[0.30] 


29.1(22.8-37.0) 


34.8 


56.3/ 45.0 


46 


IRS44/YLW16A 




2.7(2.4-3.1) 


5.3(4.9-5.7) 


0.28(0.19-0.38) 


45.6(38.3-53.9) 


49.7 


77.0/ 72.0 


66 


IRS51 




2.1(1.8-2.8) 


4.2(3.5-4.8) 


0.62(0.27-1.18) 


13.8(10.2-18.5) 


16.2 


51.7/ 48.0 



Class II 



ROXN 


NAME 




kT 


AT H 


abundance 


E.M. 


L X 


x 2 /d.o.f 


ID 








(10 22 cm" 2 ) 




(10 52 cm" 3 ) 


(10 29 erg s- 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1 


DoAr25 




6.4(4.9-8.8) 


1.1(0.96-1.2) 


0.26(0.02-0.51) 


157(141-177) 


247 


80.1/ 72.0 


11 


SR24S 




2.3(2.0-2.7) 


1.2(1.0-1.3) 


0.24(0.07-0.43) 


16.3(13.4-19.4) 


15.9 


107.7/ 75.0 


18 


WL17 




2.9(2.0-4.3) 


4.1(3.1-5.2) 


[0.30] 


8.4(5.0-14.6) 


9.6 


27.5/ 19.0 






f 


3.2(2.5-4.2) 


4.5(3.8-5.4) 


0.66(0.31-1.2) 


18.1(13.1-25.7) 


24.8 


35.4/ 25.0 


29 


IRS34 




1.8(1.1-3.7) 


2.7(1.6-4.4) 


[0.30] 


12.8(<41.5) 


11.7 


4.7/ 5.0 


30 


WL20 




2.3(1.8-3.0) 


2.0(1.6-2.4) 


[0.30] 


8.1(5.9-11.4) 


8.1 


16.2/ 17.0 






f 


3.4(2.7-4.5) 


2.0(1.7-2.3) 


0.34(0.04-0.75) 


16.4(13.3-20.2) 


20.4 


34.7/ 29.0 


33 


GY245 




2.1(1.3-5.0) 


7.7(5.2-11.0) 


[0.30] 


4.5(1.5-16.4) 


4.4 


10.9/ 9.0 


45 


VSSG25 


f 


2.0(1.6-2.5) 


1.1(0.9-1.4) 


[0.30] 


6.3(4.7-8.2) 


5.9 


29.1/ 24.0 


50 


IRS47 


f 


2.9(1.7-7.4) 


2.3(1.5-3.3) 


[0.30] 


7.5(<15.5) 


8.5 


1.2/ 2.0 


57 


GY292 




2.3(2.1-2.5) 


1.7(1.5-1.8) 


0.23(0.14-0.36) 


31.2(28.6-35.8) 


30.8 


159.3/153.0 


61 a 


IRS49 




2.0(1.9-2.2) 


1.6(1.4-1.7) 


0.18(0.08-0.29) 


27.1(23.9-30.7) 


24.1 


131.1/114.0 


62 


GY310 


f 


2.5(1.7-3.6) 


0.82(0.59-1.2) 


[0.30] 


1.8(1.4-2.4) 


1.9 


22.4/ 13.0 


63 


GY312 




2.3(1.7-3.3) 


7.6(6.0-8.7) 


[0.30] 


11.4(6.4-22.5) 


11.6 


25.0/ 17.0 


65 


GY314 




2.1(2.0-2.3) 


1.5(1.4-1.6) 


0.16(0.07-0.26) 


31.6(28.4-35.0) 


28.4 


113.2/113.0 



Class III 



ROXN 


NAME 




kT 


AT H 


abundance 


E.M. 


ix 


x 2 /d.o.f 










(10 22 cm- 2 ) 




(10 52 cm" 3 ) 


(10 29 erg s- 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


5 


GY112 




0.98(0.79-1.1) 


0.36(0.31-0.50) 


0.10(0.06-0.16) 


7.5(6.1-11.8) 


4.3 


40.9/ 37.0 


14" 


GY193 




1.27(1.19-1.74) 


0.87(0.56-1.07) 


[0.30] 


7.0(4.6-8.2) 


6.0 


29.4/ 38.0 


15 


GY194 




1.5(1.2-1.9) 


0.79(0.6-1.1) 


[0.30] 


3.2(2.6-3.9) 


2.8 


20.7/ 13.0 


25 


WL19 




2.5(1.5-4.6) 


8.62(6.0-12.8) 


[0.30] 


9.2(3.7-30.7) 


9.8 


8.4/ 10.0 


27,28 a ' 6 


ROXs20A, ROXs20B 






0.13(0.06-0.18) 


0.42(0.23-1.21) 




4.6 


59.2/ 40.0 








0.66(0.56-0.77) 






1.7(0.8-2.2) 


1.9 




32,34,35 b '° 






1.66(1.30-2.29) 






2.8(1.9-3.9) 


2.8 




WL5,WL4,WL3 




2.5(2.2-3.0) 


4.7(4.2-5.2) 


[0.30] 


60.4(47.4-76.9) 


64.2 


104.6/ 74.0 


36 


SR12A-B 






0.11(0.08-0.13) 


0.12(0.08-0.17) 




21.9 


248.6/205.0 








0.25(0.23-0.27) 






18.2(11.1-30.2) 


4.4 










0.97(0.89-1.03) 






18.6(17.0-25.2) 


11.1 










3.5(2.5-6.5) 






5.6(3.1-8.3) 


6.5 








f 




0.06(0.03-0.10) 


0.24(0.13-0.55) 




24.0 


74.0/ 60.0 








0.32(0.28-0.38) 






7.8(4.3-17.1) 


4.2 










1.22(1.09-1.35) 






15.3(5.5-22.5) 


12.0 










12.0(>3.1) 






4.9(3.0-10.4) 


7.7 




39 


GY253 




1.6(1.4-2.0) 


4.5(3.9-5.2) 


[0.30] 


24.9(17.5-34.9) 


22.0 


57.7/ 52.0 


52 


BKLT162730- 244726 




1.9(1.3-2.3) 


1.0(0.8-1.5) 


[0.30] 


6.6(5.3-11.2) 


6.1 


44.1/ 43.0 


72 


GY377 




1.8(1.1-2.7) 


1.1(0.7-1.8) 


[0.30] 


1.7(0.9-3.6) 


1.5 


7.1/ 11.0 


73 


IRS55 






1.1(1.0-1.2) 


0.31(0.20-0.50) 




62.3 


139.3/141.0 








0.64(0.60-0.68) 






46.9(28.0-67.2) 


45.7 








2.5(2.2-3.0) 






15.7(12.7-19.2) 


16.6 




76 


GY410 




1.8(1.3-2.4) 


1.5(1.1-1.9) 


[0.30] 


3.8(2.5-5.9) 


3.4 


8.9/ 23.0 


New Class III candidates 


ROXN 


NAME 




kT 


N H 


abundance 


E.M. 


ix 


x 2 /d.o.f 


ID 








(10 22 cm" 2 ) 




(10 52 cm- 3 ) 


(10 29 erg s- 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


31 


16271643-24311 




4.2(2.2-15.1) 


5.2(3.5-7.4) 


[0.30] 


6.1(3.2-14.3) 


8.3 


25.8/ 26.0 


41 


GY259 




3.4(>1.3) 


2.2(1.1-4.5) 


[0.30] 


1.9(0.8-6.3) 


2.3 


0.5/ 2.0 


44 


GY266 




2.3(1.1-6.1) 


2.8(1.6-5.0) 


[0.30] 


2.4(1.0-8.2) 


2.5 


7.5/ 6.0 


53 


GY283 




1.4(0.5-2.6) 


0.47(0.2-1.5) 


[0.30] 


0.5(0.3-3.0) 


0.5 


2.1/ 4.0 


58 


GY295 




1.9(1.2-2.8) 


0.34(0.2-0.6) 


[0.30] 


0.7(0.5-0.8) 


0.6 


12.3/ 11.0 


77 


BKLT162800-244819 




0.5(0.4-0.7) 


0.97(0.8-1.1) 


[0.30] 


8.7(4.4-17.2) 


7.8 


24.2/ 25.0 


78 


GY463 




1.6(1.2-2.1) 


2.6(2.0-3.4) 


[0.30] 


11.0(6.8-20.2) 


9.6 


16.0/ 20.0 


86 


BKLT162823-244140 




4.3(2.2-13.2) 


0.87(0.5-1.3) 


[0.30] 


4.0(2.8-6.7) 


5.5 


9.7/ 9.0 


X-ray sources 


with no optical or IR counterpart 












ROXN 


NAME 




kT 


N H 


abundance 


E.M. 


ix 


X 2 /d.o.f 


ID 








(10 22 cm- 2 ) 




(10 52 cm- 3 ) 


(10 29 erg s- 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


51 






10.5(>3.9) 


0.92(0.6-1.4) 


[0.30] 


2.5(2.0-3.8) 


4.4 


40.1/ 28.0 


64 






82.4(>20.5) 


2.2(1.7-2.9) 


[0.30] 


3.5(2.9-4.1) 


5.5 


64.5/ 43.0 


67 






3.6(>1.4) 


9.5(4.8-18.0) 


[0.30] 


6.7(2.4-46.9) 


8.6 


25.4/ 29.0 



(1)X MM- Newton source ID number, "Rho Oph X-ray sources, Newton". (2) The source name from column (15) in Table 2. (3) "f" means that the X-ray spectra 
are extracted during the flare phase. (4 — 7) The 90% confidence regions arc given between parenthesis. (5) We use WABS for absorption model in XSPEC where 
standard solar metal abundances are used. The use of rcvicsc solar abundances (Holweger 2001, Allend Prieto et al. 2001, 2001) would increase N H by 20 ~% 
(Vuong et al. 2008). (7), (8) The emission measures and X-ray luminosities were calculated from the best-fit PN parameters assuming a distance of 140 pc. Those 
from MOS1 + MOS2 parameters arc consistent within errors. 

a These sources showed small flare during the observation, but the X-ray counts in the flares arc not enough to make spectra. 

b It is difficult to obtain individual spectrum without contamination from neighbouring sources. Then we extracted X-ray spectrum from a circular region which 
contain all the sources. 

c Since the X-ray counts from the source 32 is dominant, we put the spectral parameters in the table for Class III. 
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Fig. 3. X-ray spectra of IRS44/YLW16A (Class I), EL29 (Class I), GY314 (Class II), and SR12A-B (Class III) in 
the quiescent state obtained with XMM-Newton. In each panel, the filled circles and open squares indicate PN and 
M0S1+M0S2 spectra. The solid lines show the best fit models whose spectral parameters are listed in Table 4. 



3. Comparison with Chandra source detection 

Chandra observed the p Ophiuchi cloud core F region 
with an exposure of 100 ks (Imanishi et al. 2001a), al- 
though, contrary to our XMM-Newton observations, the 
pointing was somewhat off-centered with respect to the 
peak of the CO emission of core F. Fig.^lshows the XMM- 
Newton and Chandra X-ray source positions superimposed 
on the 0.3—8 keV band image obtained by XMM-Newton. 
Chandra detected in 100 ks 87 X-ray sources in its field 
of view (17'xl7'). This the same number of sources as in 
our three times shorter XMM-Newton observations, in a 
field of view which is nearly four times as large. We shall 
return to this coincidence below. In the region covered by 
the field of view of both observatories, 47 X-ray sources are 
detected with XMM-Newton and 81 X-ray sources are de- 
tected with Chandra, while 43 X-ray sources are detected 
with both observatories. 

Fig.[S]shows XMM-Newton vs. Chandra count rates for 
all these X-ray sources. Upper limits correspond to the 
99.9989% confidence level threshold (i.e., 4.4 a level for 
Gaussian statistics) for the Chandra sources undetected 
by XMM-Newton, and for the XMM-Newton sources unde- 



tected by Chandra, using a method explained in Appendix 
A. The median of the count rate upper limits of sources 
detected with only Chandra is the dotted horizontal line 
in Fig. which indicates also the 4.4 a detection thresh- 
old for XMM-Newton. There is a good correlation be- 
tween XMM-Newton and Chandra count rates. The me- 
dian of the ratios between XMM-Newton and Chandra 
count rates for X-ray sources detected by both observato- 
ries is 6.8, which is indicated by a continuous line in Fig.0 
"the median ratio line" . Most of the data points are scat- 
tered close to this median ratio line. A few sources display 
large discrepancies between XMM-Newton and Chandra 
observations, which suggests that their X-ray luminosities 
could have changed due to time variability. We check the 
two most extreme cases, sources ROXN-33 (GY245) and 
ROXN-46 (YLW16A), which are far from the median ra- 
tio line. GY245 did not show any flare during the Chandra 
observation, but showed a flare during the XMM-Newton 
observation, while YLW16A showed an X-ray flare during 
the Chandra observation, but did not show any flare dur- 
ing the XMM-Newton observation. Hence, flaring activity 
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explains well the large apparent discrepancies observed for 
these two sources. 

There is a large difference between the number of X-ray 
sources detected with XMM-Newton and detected with 
Chandra in the region common to the field of view of both 
observatories, which we now seek to explain. The obser- 
vation time of Chandra, 100 ks, is three times longer than 
that of XMM-Newton, 33 ks. If the observation time of 
Chandra had been as short as that of XMM-Newton, the 
number of sources detected by Chandra would have been 
smaller. More precisely, if we assume that ~ 6 counts 
are needed to make a Chandra detection and that the 
count rates of the X-ray sources are constant, in a 33 ks 
Chandra observation, only those sources with count rates 
above 0.18 counts ks -1 would have been detected (see the 
dotted line in Fig. [SJ). Fig. [S] shows that most of the 23 
X-ray sources detected by Chandra only with count rates 
above 0.18 counts ks -1 are distributed well below the me- 
dian ratio line, which indicates that these sources probably 
showed X-ray flares only during the 100 ks Chandra obser- 
vation, but could disappear in a shorter, 33 ks observation. 
All in all, when the exposure time is set to about 30 ks for 
both observatories, and taking into account the differences 
in fields of view and detection thresholds, Chandra and 
XMM-Newton detect roughly the same number of sources. 

In summary, we find that the respective merits of 
XMM-Newton and Chandra, for the same type of obser- 
vations (here a typical low-mass star-forming region) are 
as follows. XMM-Newton has an effective area about three 
times larger than Chandra, so that XMM-Newton detects 
more source counts than Chandra for the same observa- 
tion time: as shown in Fig. [3 the XMM-Newton count 
rate is 6.8 times larger than the Chandra count rate. On 
the other hand, Chandra has a lower background, and a 
sharper point spread function than XMM-Newton, which 
means Chandra suffers less noise than XMM-Newton for 
point source detection. For relatively short exposures (but 
typical in a proposal request) of 30 ks, though these val- 
ues depend on position with respect to the pointing axis, 
the signal-to-noise ratios are similar for both observato- 
ries. For longer exposure times, Chandra provides a bet- 
ter detection ability than XMM-Newton, while for shorter 
exposure times, XMM-Newton provides a slightly better 
detection ability than Chandra, because the statistics with 
Chandra is photon-dominated, while for XMM-Newton it 
is background-dominated. Because it detects more pho- 
tons in a given exposure, XMM-Newton allows a more de- 
tailed spectral and time variability analysis for moderately 
bright sources (above ~ 150 counts), typical of YSOs, than 
does Chandra. 

4. IR properties of the XMM-Newton sources 

ISOCAM provides a reliable classification of YSOs in the 
p Ophiuchi cloud in terms of Class I or II sources, from 
the detection of IR excess produced by circumstellar ma- 
terial (Bontemps et al. 2001). We detect X-rays from 7/11 
Class I sources, 25/61 Class II sources, and 14/15 Class 
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Fig. 4. Comparison of source detection between XMM- 
Newton and Chandra (Imanishi et al. 2001a). The back- 
ground image shows the XMM-Newton image in the 0.3—8 
keV energy range smoothed with a 6"-FWHM Gaussian. 
The rectangle shows the field of view of Chandra. The tri- 
angles, circles, and crosses show the positions of the X-ray 
sources detected with only XMM-Newton, only Chandra, 
and both observatories, respectively. 

Ill sources in the overlapping region between the field of 
view of XMM-Newton and the ISOCAM survey. 7 We find 
a high detection rate for Class I sources (64 %), which is 
higher than that of Class II sources (48 %). The XMM- 
Newton detection rate of Class I sources is comparable to 
that obtained from the Chandra observation, but the de- 
tection rate for Class II sources is lower by more than 70 
% than that obtained by Chandra (Imanishi et al. 2001a). 
The Class III sample is biased because these sources were 
mostly identified by X-rays, which explains why their high 
detection rate (93 %). 

Fig. El is a scatter plot of stellar luminosities (L*) vs. 
extinction (Ay), for Class II (filled circles) and Class III 
sources (open circles) from Bontemps et al. (2001) results. 
Sources having XMM-Newton X-ray counterparts are in- 
dicated by crosses. Most of the X-ray non-detected Class 
II sources are distributed in the lower-right region (low 
and high Ay) of Fig.El Assuming that there is an Lx / L+ 
correlation (see below, Sect. 6), combined with the high 
extinction, the X-ray non-detected Class II sources in this 
region could simply lie below the detection threshold of 
XMM-Newton (see Grosso et al. 2000) which, as discussed 
above, is higher than for Chandra. 

However, some X-ray non-detected Class II sources, 
like SR24N , WSB37, and GY3, are not located in the 
lower-right region of Fig. [S] This could be explained in two 

7 Note that the prototype of Class sources, VLA1623 
(Andre et al. 1993) is located in the dense core A, which is 
not covered by our observation. 
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Fig. 5. Comparison between XMM-Newton/EPIG and 
Chandra/ 'ACIS-I count rates for X-ray sources in Chandra 
and XMM-Newton overlapping ficld-of-views. The contin- 
uous line shows the median of the count rate ratio for 
sources detected by both instruments. The dashed line 
is the median of the count rate upper limit for Chandra 
sources which are not detected during the 33 ks XMM- 
Newton observation. The dotted line works the Chandra 
detection limit (~ 6 counts) for a 33 ks observation. 



ways. First, their individual Lx/L+ ratios may be lower 
than typical ones. Second, their extinction may be higher 
than given by the statistical extinctions estimated from 
the NIR data (Bontemps et al. 2001). This could happen if 
the circumstellar disk of these Class II sources is seen close 
to edge-on: in such a case, the central star is not visible, 
and we just observe its scattered light so that a statisti- 
cal value of Ay underestimates the real value. Conversely, 
two luminous Class III sources with very high extinction, 
WL5 (ROXN-32) and WL19 (ROXN-25), are located in 
the upper-right corner of Fig. H3 Their X-ray spectral fits 
provide a column density consistent with high absorption 
(see Table 4). These two Class III sources are probably 
located behind the dense core F of cloud. A similar con- 
clusion, in the much wider field of view of ROSAT, was 
reached for many Class III sources by Grosso et al. (2000). 
The correlation between X-ray absorption column density 
and optical extinction was studied in detail by Vuong et al. 
(2003), who compare the gas and dust properties of the 
dense interstellar matter in nearby star-forming regions 
using X-rays from Class III sources. 

Considering that ISOCAM observations of Bontemps 
et al. (2001) were sensitive enough to detect IR excess, if 
any, from the new X-ray sources found in our observation, 
they could be Class III candidates. There are 37 XMM- 
Newton X-ray sources without IR classification in Table 2, 
18 of which have 2MASS counterparts. With the exception 
of ROXN-10 (identified with the bona fide brown dwarf 
GY 141, see below, §5), we check now whether the near- 
IR photometry of these X-ray sources is consistent with 
that of Class III sources. 



Fig. [7| shows the color-color diagram of the XMM- 
Newton X-ray sources identified with the 2MASS near- 
IR sources. Class I, II, and III sources and brown dwarfs 
are indicated by filled diamonds, filled circles, open cir- 
cles, and open diamonds, respectively; and sources with- 
out IR classification are indicated by asterisks labeled with 
the ROXN numbers of Table 2. As ROXN-31 is too faint 
to have measurable J and H magnitudes, it is removed 
from our sample, and 16 unclassified sources are plotted 
in Fig.0 For ROXN-41, -44, and -49, only lower limits are 
available for the J — H color. ROXN-87, corresponding to 
ROSAT source ROXR-F37, was identified as a foreground 
F2 V star, HD148352 (Grosso et al. 2000), which is con- 
sistent with its position on the locus of the intrinsic colors 
(Bessel & Brett 1988). This is also consistent with the 
fact that its X-ray hardness ratio is -1.0 for MOS1 and - 
1.0 for MOS2 from Table 1, which means there is no hard 
X-ray emission above 2 keV. ROXN-47 is also located at 
the position of a late-type M dwarf without extinction in 
Fig. its hardness ratio is relatively low, —0.54, —1.00, 
and -0.76 for MOS1, MOS2, and PN, respectively. We 
thus consider this source also as a foreground star 

Fig. |S| shows the color-magnitude diagrams of H vs. 
J — H in the left panel, and K s vs. H — K s in the right 
panel for XMM-Newton sources having 2MASS counter- 
parts. (Note that ROXN-31 appears only in the K s vs. 
H — K s diagram since this source is detected only in the 
K s band.) We plot for comparison the 1 Myr isochrone 
(Baraffe et al. 1998) and reddening vectors (Cohen et al. 
1981). The positions of the remaining 15 sources without 
IR classification are well mixed with those of well-known 
Class II and III sources. Therefore, we propose these 15 
X-ray sources as new Class III candidates (labeled 'nlll' in 
Table 2), therefore doubling the present number of Class 
III sources in this area. Spectroscopic follow-up is now 
needed to determine the effective temperature of these 
objects, and to put them in an H-R diagram to confirm 
their pre-main sequence status. We estimate stellar lumi- 
nosities and extinctions of these new Class III candidates 
using J and H band 2MASS magnitudes (see formulas in 
Bontemps et al. 2001), and plot them in Fig. The new 
Class III candidates are also mixed with the Class II and 
III sources in Fig. [5] 

In summary, we find a high XMM-Newton detection 
rate for Class I and Class III sources, consistent with the 
Chandra results, and a lower detection rate for Class II 
sources. For those sources, the difference is probably due 
to a combination of high extinction (interstellar + circum- 
stellar) and of lower detection efficiency of XMM-Newton 
compared with Chandra. We however detect 15 previously 
unknown X-ray sources, which we propose as new Class 
III candidates, pending spectroscopic follow-ups to con- 
firm their nature. Our XMM-Newton observations thus 
allow for a significant improvement of the YSO census in 
the p Oph cloud core F region (15 new YSOs in addition 
to a total of 87 previously known from X-ray/IR observa- 
tions, and a potential doubling of the number of Class III 
sources) . 
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Fig. 8. Color-magnitude diagram of the XMM-Newton X-ray sources with 2MASS counterparts. The solid line shows 
the 1 Myr isochrone for solar metalicity low mass stars from 0.02 to 1.00 Mq (Baraffe et al. 1998), and 10 mag 
reddening vectors (Cohen et al. 1981). Symbols are the same as in Fig. Arrows indicate lower limits on J — H and 
H — K a when J and H magnitudes are not available, respectively. 
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Fig. 6. Stellar luminosity vs. extinction for Class II (filled 
circles) and III (open circles) sources in the ISOCAM sur- 
vey (Bontemps et al. 2001) and XMM-Newton overlapping 
area. The sources which have XMM-Newton counterparts 
are indicated by crosses. The new Class III candidates are 
indicated by asterisks and are labeled with ROXN num- 
bers of Table 1 and 2. 



5. X-ray detection of young bona fide brown 
dwarfs 

Since ROSAT observations, it is known that young brown 
dwarfs also emit X-rays (Neuhauser & Comeron 1998; 
Neuhauser et al. 1999; Imanishi et al. 2001a; Imanishi et 
al. 2001b; Preibisch & Zinnecker 2001; Mokler & Stelzer 
2002; Feigelson et al. 2002b; Tsuboi el al. 2003). We thus 
look for X-rays from young bona fide brown dwarfs in our 
observation, i.e., from young objects with substellar status 
confirmed by spectroscopy. There are in our field-of-view 
only three bona fide brown dwarfs 8 : GY202, GY141, and 
GY310 (Martin et al. 1999; Cushing et al. 2000; Wilking 
et al. 1999). Bontemps et al. (2001) classified GY310 as 
Class II, and recently Mohanty et al. (2004) detected it 
from the ground with Subaru at 8.6 and 11.7/im, confirm- 
ing the presence of significant mid-infrared excess arising 
from an optically thick, flared dusty disk. 

A low S/N ratio X-ray detection of GY202 was re- 
ported by Neuhauser et al. (1999) from the ROSAT/PSPC 
pointing observation of Casanova et al. (1995). We find 
that the identification of this ROSAT/PSPC source with 
GY202, located 19" away, is dubious because the closest 
counterpart is in fact WL1, an embedded (Ay ~ 21 mag) 

8 We note in Fig.^that ROXN-82 is located well below the 
substellar limit of ~ 0.08 Mq and could be a brown dwarf 
candidate, but it was found to be a very low-mass M star by a 
near-IR spectroscopical follow-up (Grosso et al. 2004, in prepa- 
ration) . 
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Fig. 7. Color-color diagram of the XMM-Newton sources 
which have 2MASS near-IR counterparts. The arrow 
shows an extinction of 10 mag (Cohen et al. 1981). The 
intrinsic colors of giants (dotted line) and AO— M6 dwarfs 
from Bessel k Brett (1988), adapted for the 2MASS pho- 
tometric system using the 2MASS color transformation 
(Carpenter 2001; Cutri et al. 2003), are plotted for com- 
parison. The dashed line shows the locus of the intrinsic 
color classical T Tauri stars (Meyer et al. 1997). Arrows 
indicate lower limits on J — H when J magnitude is not 
available. 



Class II source (Bontemps et al. 2001), located only 
9" away. Moreover, in spite of the fact that we do not de- 
tect WL1, its X-ray emission is confirmed by the Chandra 
observation (source 13 of Imanishi et al. 2001a) with a 
luminosity roughly consistent with the ROSAT/PSPC es- 
timate, whereas GY202 is not detected either by Chandra 
(Imanishi et al. 2001b; Imanishi et al. 2003) and XMM- 
Newton (this work). 

A very weak X-ray emission has been reported from 
GY141 with Chandra (Imanishi et al. 2001b; source BF-S2 
in Imanishi et al. 2003): only ^8 X-ray photons were col- 
lected during the 100 ks exposure. According to Fig. lathis 
low count rate is well below the sensitivity of our XMM- 
Newton observation. However we detect GY141 (ROXN- 
10) with a count rate of 7.7ctsks~ 1 , i.e., at a level ~ 
90 times higher than during the Chandra observation. 
Although we detect this brown dwarf during a phase 
of intense X-ray activity, we do not have enough statis- 
tics to investigate further this high X-ray state (see its 
background-subtracted light curve in Fig. [[JJ. 

X-ray emission was also reported from GY310 with 
Chandra (Imanishi et al. 2001a; Imanishi et al. 2001b; 
Imanishi et al. 2003). The Chandra X-ray light curve 
shows no clear flare, but exhibits aperiodic variability 



by a factor 2 within the 100 ks exposure, around ~ 2 x 
10 _3 ctss _1 (Imanishi et al. 2001b). During our obser- 
vation, GY310 (ROXN-62) displayed an X-ray flare (see 
Fig-El- To our knowledge this is the first X-ray flare from a 
young bona fide brown dwarf with enough counts to derive 
its X-ray spectrum due to XMM-Newton's large effec- 
tive area. The observed count rate with MOS1+MOS2 
increased from 0.007ctss _1 to 0.04ctss _1 . The quies- 
cent level gives an equivalent Chandra count rate of 
2 x 0.007/6.8 - 2 x lO^ctss -1 , consistent with the 
low value previously observed by Chandra. We obtain 
enough X-ray counts from GY310 during the flare to make 
spectral analysis presented in Sect. 12.31 (see Fig. ^| and 
Table 4). The derived column density is identical to the 
one found by Imanishi et al. (2001b). We find a some- 
what higher plasma temperature in our flare observation 
(2.5 keV, 1.7-3.6 keV), compared to the quiescent Chandra 
value (1.7 keV, 0.9-2.2keV), a fairly general characteristic 
of stellar X-ray flares. 

Similar changes in the level of X-ray activity, X-ray 
flares, and high plasma temperatures, are ubiquitous in 
low-mass protostars and T Tauri stars. This suggests that 
X-ray emission from brown dwarfs in their early phase of 
evolution are produced basically by the same solar-like, 
magnetic activity mechanism at work in low-mass proto- 
stars and T Tauri stars, and more generally in late-type 
stars, as also noticed in previous studies (Imanishi et al. 
2001b, Feigelson et al. 2002b). 

6. X-ray properties of Young Stellar Objects in 
the p Ophiuchi cloud 

To investigate a possible correlation between the X-ray 
luminosity Lx and the stellar luminosity of YSOs, we 
plot Lx vs. for Class I, II and III sources, new Class III 
candidates, and brown dwarfs in Fig.^] We calculate ab- 
sorption corrected Lx in the 0.5—10 keV band using the 
best fit model for the X-ray spectra, which are listed in 
Table 4. In Fig.f^ we plot only the ROXN sources which 
are bright enough to determine their spectral parameters, 
except for ROXs20A (ROXN-27), ROXs20B (ROXN-28), 
WL5 (ROXN-32), WL4 (ROXN-34), and WL3 (ROXN- 
35), for which the individual X-ray spectra cannot be re- 
solved. We use L» from Bontemps et al. (2001) for Class 
I, II, and III sources and brown dwarfs, and use our es- 
timation of L* for new Class III candidates presented in 
section 0] For Class I sources, as the stellar luminosity is 
unknown because the central star is invisible due to the 
remnant dust envelopes, we use the bolometric luminosity 
as an upper limit to the stellar luminosity. The correlation 
index between Lx and L* of Class II and III sources in 
the quiescent state is then —0.08, which indicates a weak 
correlation. Indeed the Lx / L* ratios show a large spread, 
from 4.5 x 10~ 6 to 1.6 x 10~ 2 . New Class III candidates 
have similar Lx and L* properties as other YSOs, which 
is consistent with the assumption of their YSO nature. 

To compare the characteristics of the X-ray spectra 
of Class I, II, and III sources, we show in Fig. ffJl the 
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Fig. 9. X-ray background subtracted light curves of the 
young bona fide brown dwarfs, GY310 and GY141, ob- 
tained with XMM-Newton. High background time inter- 
vals were suppressed (holes in the light curve). We use an 
adaptive binning to keep a constant count number in each 
bin (see Table 3). 



scatter plot of their X-ray determined absorption column 
density, TVh, vs. plasma temperature, kT. For ROXN-36 
(SR12A-B) and ROXN-73 (IRS55), we calculated the av- 
erage temperature from the multiple components weighted 
by their emission measures (see above, §2.3). An interest- 
ing property emerges from the JVh vs. kT diagram: only a 
few data points are seen in the upper-left region and the 
lower-right region. The upper-left region is where the X- 
ray sources have a low plasma temperature and large ab- 
sorption: since the soft X-rays from low temperature plas- 
mas are easily absorbed, it is natural that X-rays should 
be hard to detect from these sources. On the other hand, 
the lower-right region is where the X-ray sources have a 
high plasma temperature and small absorption, and are 
thus easy to detect, but this region contains surprisingly 
few data points. This could be explained by an intrin- 
sic effect if both the maximum temperature of YSOs and 
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Fig. 10. X-ray spectra of the young bona fide brown 
dwarf GY310 during the flare state obtained with XMM- 
Newton. The filled circles and open squares indicate PN 
and MOS1+MOS2 spectra. The solid lines show the best 
fit models whose spectral parameters are listed in Table 
4. 

their circumstellar material causing absorption of X-rays 
decrease in the course of their evolution. Class I sources, 
which arc in the early stage of the evolution, arc indeed 
located in the upper-right region of this diagram. Class II 
and III sources, which are in a later stage of the evolution, 
are located in the region of lower temperature and smaller 
absorption than Class I sources, which confirms the results 
obtained by Chandra (Imanishi ct al. 2001a). There is also 
a tendency for the temperature and absorption of Class III 
sources to be lower than those of Class II. We conclude 
that, although our statistics are still limited, an evolution- 
ary effect seems to be present from the Class I stage (high 
X-ray temperature, high extinction) to the Class III stage 
(lower X-ray temperature, lower extinction). 

7. Summary and conclusions 

The main results of our XMM-Newton 33 ks, 30' diameter 
observation of the p Ophiuchi cloud core F (i.e., the main 
molecular core) are as follows : 

1. We detect 87 X-ray sources; 62 sources are already 
known from previous X-ray observations, hence we find 
25 new X-ray sources. 

2. In the region covered by the field-of-view of both 
XMM-Newton and Chandra observations, 4 X-ray 
sources are detected only with XMM-Newton and 38 
X-ray sources are detected only with Chandra, while 
43 X-ray sources are detected with both observatories. 
We have shown that these differences are explained by 
the difference in exposure time and sensitivity to point 
source detection, along with intrinsic time variability 
of X-ray sources (flares). We find that Chandra and 
XMM-Newton are equally sensitive for point source de- 
tection in star forming-regions for exposures of <~ 30 
ks. 
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Fig. 11. X-ray luminosity, Lx, vs. stellar luminosity, L*, 
plot for YSOs in the p Ophiuchi cloud core F. The open 
squares, stars, triangles, and diamonds indicate Class I, 
II, and III sources and brown dwarfs in the flare state, re- 
spectively, and the filled squares, stars, triangles, and cir- 
cles indicate Class I, II, and III sources and new Class III 
candidates in the quiescent state, respectively. For Class 
I sources, the bolometric luminosities are plotted as the 
upper limits of stellar luminosities. 
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Fig. 12. Nn vs. kT plot of the X-ray sources listed in 
Table 4. The symbols are the same as in Fig. ^2 For 
ROXN-27-28, ROXN-36, and ROXN-73 which need multi- 
temperature model for the spectral fitting, the average 
temperatures weighted by the emission measures are plot- 
ted. The errors show 90 % confidence region. 



3. We detect X-ray emission from 7 Class I sources, 26 
Class II sources, and 17 Class III sources. For Class I 
sources, our detection rate is high (64 %), consistent 
with that already observed by Chandra (Imanishi et 
al. 2001a). 



For Class II sources, our detection rate is much lower 
(48 % vs. 70 %), which can be explained by a combi- 
nation of instrumental reasons (better sensitivity for 
point source detection by Chandra) and intrinsic rea- 
sons (better probability to detect a source via X-ray 
flares in a longer exposure). 

We find a high detection rate for Class III sources, 
which is not surprising in view of the fact that most 
Class III sources were discovered in X-rays. But we 
propose 15 previously unknown X-ray sources as new 
class III candidates, in addition to the already known 
87 YSOs in the p Oph core region (Bontemps et al. 
2001). Pending spectroscopic confirmation, this would 
double the number of Class III sources known in this 
area. 

6. We detect X-ray emission from two young bona fide 
brown dwarfs, GY310 and GY141. GY310 showed 
an X-ray flare during our observation, and GY141 
appeared brighter by nearly two orders of magni- 
tude than in the Chandra observation (Imanishi et al. 
2001a). Their X-ray properties suggest that the X-ray 
emission from young brown dwarfs is the same as in 
low-mass protostars and T Tauri stars, i.e., resulting 
from solar-like magnetic activity, as proposed by pre- 
vious studies (Imanishi et al. 2001b, Feigelson et al. 
2002b). 

7. We are able to extract X-ray light curves and spectra 
from a dozen sources, and find some of them showed 
X-ray flares. Altogether, from spectral fitting we con- 
firm that there is an evolutionary trend, from Class I 
sources showing higher temperature and larger absorp- 
tion, to Class II and III sources showing lower temper- 
ature and smaller extinction, as previously reported by 
the Chandra observation (Imanishi et al. 2001a). 
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exposure map with a spatial resolution of 2" for X-ray 
photons of 1.5 keV energy, and normalize it by its maxi- 
mum value located on-axis. The resulting bidimensionnal 
array, £, represents the spatial effective area variations, 
including CCD gaps, bad pixels, and vignetting, relative 
to the on-axis value. From £ we define an exposure mask, 
M, having M(i)=0 if£(i)=0, and M(i)=l if£(i)^0. This 
exposure mask shows only CCD gaps. 

We estimate for each position of undetected Chandra 
source the local average background. We extract the event 
number, iVbgd> inside the 15"-radius disk, 2?b g d, of geomet- 
rical area centered on the Chandra source position. 
This extraction region contains n^ gd sky pixels of £ , and 
n bgd pixels of M where M is equal to 1, i.e. which 
are not CCD gaps. The net pixel area of Pbgd is : 

^nct _ Ageo *M I £ 
^bgd — ^bgd n bgd / n bgd- 



(A.l) 



The average local background per sky pixel is obtained 
by: 



<B>=N hgd /A^ d . 



(A.2) 



This 15"-radius size is a compromise for crowded re- 
gion. With this local method, local background variations 
as PSF wings of neighbouring X-ray source, or trails due to 
out-of-time events of bright sources, are automatically in- 
cluded in our background estimate. This would not be the 
case if we had used the spline-smoothed background built 
by esplinemap, because it uses a 'Swiss-cheese' method 
which subtracts the source PSF from the X-ray image. 

Then we compute the PSF radius encircling / cc = 50% 
of the PSF energy, i?(/ ee ), using the radial average of the 
calibration images in the CCF. We extract the event num- 
ber, A^src+bgdC/co), within the _R(/ oc )-radius disk, X> src , of 
geometrical area centered on the Chandra source po- 
sition. These counts are the sum of background and source 
counts, A r src (/ ee ) and -/Vb g d(/ee), respectively. The extrac- 
tion region Aff° contains nf rc sky pixels of £, and n*^ 
sky pixels of M where M is equal to 1. The net pixel area 
of X> src is : 



\ net 



A goo *M 
^sre ,f, src 



£ 

sre* 



(A.3) 



We compute for iV src+ b g d(/cc) its 99.9989% confidence 
level (corresponding to 4.4(7 in Gaussian statistics, i.e. 
the threshold of our detection algorithm) upper limit, 
N^ c+h d (f ee ), using the approximate formula (10) of 
Gehrels (1986). Assuming that the background is constant 
on V src , the upper limit on the source count in T> src is : 



ATr'cC/ee) 



{<c + bgd(/ee)-<S>^ r C c t } 



Appendix A: Method to compute the X-ray count 
rate upper limit 

A.l. XMM-Newton data 

For each instrument, MOS1, MOS2, and PN, we com- 
pute with the SAS command eexpmap the corresponding 



(A.4) 



The right hand side multiplicative factor corrects in aver- 
age from loss of sensitivity due to bad pixels, CCD gaps, 
and vignetting. For Chandra sources which are over CCD 
gaps or bad pixels of XMM-Newton (a minority in our 
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sample), a more accurate correction would involve a con- 
volution of the PSF shape with the exposure mask. 
The upper limit on the total source count is hence : 

N*=N£ c (f ee )/f ee . (A.5) 

Finally, the count rate upper limit is obtained by dividing 
N^rc by the livetime of the corresponding CCD. 

We compute the count rate upper limit for the sum 
of the MOS1, MOS2, and PN data, the so-called EPIC 
instrument, CR^. C EPIC , using this set of straightfor- 
ward formulae using the previous notation, where i G 
[MOSl,MOS2,PN] : 

N src+bgd,EPIC 

(/ ee ) = AWbgd,i(/eo) , (A.6) 



-^src,EPIc(/ee) — j^c+bgcut/cc) ~ E < B >' A "c,i j 

x ^ sre (A.7) 

•^src.EPIC = ^src,EPIc(/cc) / fee , (A-8) 

CKUpic = K'cEPic / E LIVETIME, . (A.9) 

i 

A.2. Chandra data 

We get the level 2 data of the Core F observation of the p 
Ophiuchi dark cloud (sequence number 200060, observa- 
tion ID 635) from the Chandra data archive 9 . We compute 
the exposure map with a resolution of 4" for X-ray photons 
of 1.5 keV energy, and normalize it by its maximum value 
(on- axis), from which we derive an exposure mask. We 
proceed as for XMM-Newton data, using for the PSF ra- 
dius encircling 50% of the PSF energy the formula given by 
Feigelson et al. (2002b) : i?(50%) = 0.43 - 0.1 6 + 0.05 6> 2 , 
where 9 is the off-axis angle in arcmin. 

We note that for Chandra data, thanks to the large 
satellite wobbling, bad pixels and CCD gaps are smooth 
in the exposure map, they do not produce holes with sharp 
edges, hence with the above notation we have n* M /n £ = 
1. By combining (l)-(3), we find a formula identical to the 
formula (6) of Feigelson et al. (2002b) when replacing the 
90% upper limit count by the observed count. 




N n (10 22 cm 2 ), XMM-Newton 

Fig. B.l. Comparison of Ah determinations from spec- 
tral fitting between Chandra and XMM-Newton. We plot 
for comparison the diagonal with a slope equal to 1. The 
symbols are the same as in Fig. 10. The errors show 90 % 
confidence region. Chandra and XMM-Newton estimates 
of Ah are consistent taking account the error bars. 

IB. If) . A best fit value of the ratio of the iVn values from 
Chandra and XMM-Newton obtained from a linear fit is 
0.96 ( 0.93—0.99, 90 % confidence region), indicating that 
the both determinations of Ah values are consistent better 
than 10%. 

Uncertainties of the X-ray derived Ah values were dis- 
cussed for both Chandra and XMM-Newton in Vuong et al. 
(2003). They showed that the use of the recently revised 
solar abundances (Holweger 2001, Allend Prieto et al. 
2001, 2002) in X-ray spectral fitting increases N-^ values 
by ^20 %. Both in the X-ray spectral fitting by Chandra 
(Imanishi et al. 2001a) and by XMM-Newton (this work, 
see Sect. l2.3|) . the Ah values were determined using WABS 
absorption model in XSPEC where solar standard metal 
abundances are used. If the abundance effect shown by 
Vuong et al. (2003) is taken into account, all the Ah val- 
ues derived by Chandra and XMM-Newton would increase 
and all the data points in Fig. IB. II would move towards 
upper-right direction by ~20 % along the diagonal. 



Appendix B: Comparison between Chandra and 
XMM-Newton determination of Ah 

To check the reliability of the X-ray derived values of 
Ah, we make a plot of the Ah values determined by 
Chandra vs. those determined by XMM-Newton for the 
X-ray sources with spectral fitting data (see Table. 4) in 
the overlapping field-of- views of the two satellites (see Fig. 

9 http://asc.harvard.edu/cgi-gen/cda/retrieve5.pl . 



